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Abstract. This paper surveys various aspects of the hydrodynamic formula- 
tion of the nonlinear Schrodinger equation obtained via the Madelung trans- 
form in connexion to models of quantum hydrodynamics and to compressible 
fluids of the Kortcwcg type. 



1. Introduction 

In his seminal work [371 i^'^'^ ^1^° [33]): E. Madelung introduced the so-called 
Madelung transform in order to relate the linear Schrodinger equation to a hydro- 
dynamic type system. This system takes (slightly) different forms according to the 
context: linear or nonlinear equation with various nonlinearities. For the semi- 
classical nonlinear Schrodinger equation (shortened in NLS in what follows): 

^e^t^P' + ^-^Ar = f{\r^ ^^ 
the Madelung transform amounts to setting 

so as to get the following system for p and v = V^: 



(1.1) 



dtv + v-Vv + Vf{p) = — V 



/A(Vp) 



^ \ VP J 

^ dtp + div (pv) = 0. 

This system is referred to as the hydrodynamic form of NLS because of its sim- 
ilarity with the compressible Euler equation (which corresponds to e—O). The 
additional term on the right-hand side is the so-called quantum pressure. 

Madelung transform is crucial to investigate qualitative properties of the non- 
linear Schrodinger equation with nonzero boundary conditions at infinity whenever 
the solution is not expected to vanish "too often". Of particular interest is the 
so-called Gross-Pitaevskii equation 

+ yA^^ = dVl' - 1)V^ 

which corresponds to /(r) = r — 1 or, more generally, the case where / is a smooth 
function vanishing at some and such that /'(j^o) < 0. This covers in particular 
the "cubic-quintic" NLS {f{r) = —ol\ + a-^r — a^r'^ with ai, as, as > 0). 

Using Madelung transform is also rather popular to study the semi-classical limit 
of NLS. This requires p to be nonvanishing, though. We shall see below to what 
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extent the presence of vacuum (that is of points where p vanishes) is merely a 
technical problem due to the approach related to Madelung transform. 

The hydrodynamic form of NLS may be seen as a particular case of the system 
of quantum fluids (QHD) with a suitable choice of the pressure law, which, in 
turn, enters in the class of Korteweg (or capillary) fluids. We aim here at further 
investigating the link between these three a priori disjoint domains, through various 
uses of the Madelung transform. In passing, we will also review some more or less 
known facts pertaining to the theories of quantum fluids and Korteweg fluids. The 
most original part of this work is a simple proof of a recent result of AntoncUi and 
Marcati [31I3|, on the global existence of weak solutions to a quantum fluids system. 

1.1. Organization of the paper. The paper is organized as follows. The second 
section is a review of the connexions of the Madelung transform with the semi- 
classical limit of the nonlinear Schrodinger equation. In passing we briefly present 
the state-of-the-art for the Cauchy problem for the Gross-Pitacvskii equation. In 
Section [3l we solve the quantum hydrodynamical system (jl.ip by a direct method 
based on the use of an extended formulation, and explain how it may be adapted 
to tackle general Korteweg fluids. In Section 21 we review the use of the Madelung 
transform and of the hydrodynamical form of the Gross-Pitaevskii equation to 
study the existence and properties of its traveling wave solutions and of its tran- 
sonic limit, both in the steady and unsteady cases. We give in Section [5] a simple 
proof of the aforementioned result of Antonelli and Marcati ([Sill]). Lastly we list 
in Appendix the basic conservation laws for the Schrodinger, the QHD and the 
compressible Euler equations, and show that some of these laws naturally carry 
over to general Korteweg fluids. 

1.2. Notations. 

• Wc denote by (1 ^ p ^ oo) the standard norm of the Lebesgue spaces 
LP(R'*). 

• The standard H^'iR'^) Sobolev norm will be denoted 11 • lie. 

• Wc use the Fourier multiplier notation: f{D)u is deflned as F{f{D)u){S) ~ 
/(^)m(^), where J- and ^ stand for the Fourier transform. 

• The operator A = (1 — A)^/^ is equivalcntly defined using the Fourier 
multiplier notation to be A = (1 + |_Dp)^/^. 

• The partial derivatives will be denoted with a subscript, e.g. ut, u^, ... or 
dtu, dxU, ... or even dj (to designate dx^). 

• C will denote various nonnegative absolute constants, the meaning of which 
will be clear from the context. 

2. Madelung transform and the semi-classical limit of NLS 
As already mentioned in the Introduction, we consider the equation 

(2.1) isdt^ + = ; V^nO,x) = v/^e*'*"^^)/^ 

Here we have in mind the limit £ — ?> 0. The space variable x belongs to in 
this paragraph. The periodic case x G T'' being of particular interest for numerical 
simulations however (see e.g. [SI US] and references therein), we will state some 
analogous results in that case, which turn out to be a little simpler than in the case 
X G W^. In this section, we shall focus on two types of nonlinearity: 
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• Cubic, defocusing nonlinearitj0: 

f{\r\') ^^ = iviv. 

The Hamiltonian associated to ()2.1|) then reads 

m = \\e^r\\h + u'\\l^- 

It is weh defined on the Sobolev space in dimension d ^ 4. 

• Gross-Pitaevskii equation: 

f{\i,^^)r = {\r\'-i) 

In this case, the natural Hamihonian associated to (|2.ip is 

nhp m = \\svi^'\\i. + \\m'-A\l.- 

In both cases, the Hamiltonian defines an energy space, in which existence and 
uniqueness for the Cauchy problem p.ip have been established. 

The case of a defocusing cubic nonlinearity is now well understood. In dimension 
ti ^ 4 the corresponding NLS equation is globally well-posed in H^{R'^), and the 
additional iJ'*(M'') regularity (s > 1) is propagated (see the textbooks [T8l [36l |44] 
and the references therein). 

The situation is more complicated for the Gross-Pitaevskii equation, where the 
finite energy solutions ip^ cannot be expected to be in L^(R'^), since IV'^P ~ 1 G 
L'^{M.'^). As noticed in [49l [50], and extended in [24], a convenient space to study 
the Gross-Pitaevskii equation is the Zhidkov space: 

X''{R'^) = {i' e L°°{R'^) ; VV' e H''~\R'^)} with s > d/2. 

Remark 2.1. In the case x G T"^', the spaces H-'{T^) and X-'{T') (with obvious 
definitions) are the same. 

In the case a; € M'', as a consequence of the Hardy-Littlewood-Sobolev inequality 
(see e.g. [35J Th. 4.5.9] or Lemma 7]), one may show that if d ^ 2 and 
^ e I?'(E'^) is such that e LP{R'^) for some p e]l,d[, then there exists a 
constant 7 such that ip ~ 1 & L'^(R'^), with 1/p = 1/q + 1/d. Morally, 7 is the 
limit of ip at infinity. If d ^ 3 then we can take p = 2, so every function in X^(M'^) 
satisfies the above property; for d = 2, the above assumption requires a little more 
decay on Wip than general functions in X'*(M^). 

The well-posedness issue in the natural energy space 

associated to the Gross-Pitaevskii equation has been investigated only very recently 
by C. Gallo in [25] and P. Gerard in [28l [29], in dimension d ^ 4. 

Let us emphasize that in the M'' case, the energy space is no longer a linear space 
(contrary to the case of zero boundary condition a infinity where it is iJ^(R")), 
hence solving the Gross-Pitaevskii equation in iJ(R'^) is more complicated. How- 
ever, if d = 3, 4, one may show that E{M.'^) coincides with 

= c(l + v),ce S\tp e H\R'^),2 Re^/j + £ L^{R'^)}, 

which allows to endow it with a structure of a metric space. The case d = 2 is 
slightly more technical. 



^We could consider more general defocusing nonlinearities such as / (|V''^P) V''^ = li'^l^" ''P'^ t 
(T S N. We consider the exactly cubic case for the simplicity of the exposition only. 
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2.1. Some issues related to the use of the Madelung transform. In the 

semi-classical context, Madelung transform consists in seeking 

(2.2) ^P'{t, x) = Vp(MOe"^^*'''-'^^ 

for some p ^ and real-valued function cj). Plugging (j2.2p into ()2.ip and separating 
real and imaginary values yields: 

j ^(^dt^+^\Vc^\'+fip)^ ^'-^A^ ; 0|t=o = 0o, 
(2-3) < ^ 

Two comments are in order at this stage: the first equation shows that (f> depends on 
e and the second equation shows that so does p in general. We shall underscore this 
fact by using the notation (0'^,^^). Second, the equation for (j)'^ can be simplified, 
provided that p'^ has no zero. Introducing the velocity v'^ = ^(jf , (|2.3p yields the 
system of quantum hydrodynamics (jl.ip presented in the introduction. 

To study the limit £ ^ (the Euler limit), it is natural to consider the following 
compressible Euler equation: 




ftw -t- w • Vw + V/ (p) = ; W|t=o = V(/>o, 
dtp + AiY{pv) ; p\t=o^PQ- 



Note that the solution to p.4p must not be expected to remain smooth for all time, 
even if the initial data are smooth. In [38] (see also [47]), it is shown that compactly 
supported initial data lead to the formation of singularities in finite time. In [42] , the 
author constructs a solution developing singularities in finite time, in the absence 
of vacuum. In [2], it is shown that for rotationally invariant two-dimensional data 
that are perturbation of size e of a rest state, blow-up occurs at time Tg ^ t je^ . 

Remark 2.2. In the framework of this section, we have /' = 1 > 0. As pointed 
out above, the case /(r) = r'^, cr e N, could be considered as well. On the other 
hand, if /' < (corresponding to the semiclassical limit for a focusing nonlinearity), 
p.4p becomes an elliptic system which may be solved locally-in-timc for analytic 
data (see [27l |46]). At the same time, in the case d —\ and /' = — 1. it has been 
shown in [40] that there are smooth initial data for which the Cauchy problem (|2.4|) 
has no solution. In short, working with analytic data in this context is not only 
convenient, it is mandatory. 

Let us emphasize that the hydrodynamical formulations for the cubic NLS and 
the Gross-Pitaevskii equations are exactly the same as in both cases we have 
V/(p) = Vp. From this point of view, studying either of the equations is mainly a 
matter of boundary conditions at infinity: is the appropriate space for the 

cubic NLS equation whereas Ar''(IR'') is adapted to the Gross-Pitaevskii equation. 

In the sequel, denotes cither i/''(]R'') or A'''(]R''). In addition, we set 



s>d/2 
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Theorem 2.3. Let po,(j)o G C°°(M'') with y/p„,'^(f)o e Z'' for some s > d/2 + 1. 
There exists a unique maximal solution {v,p) e C([0, Tmax); ■Z^'*) to (j2.4p . In addi- 
tion, Tijiax independent of s > d/2 + 1 and 



Jo 

Finally, if V0o ^.i^d po are smooth, nonzero and compactly supported then Tmax is 
finite. 

We investigate the following natural questions: 

Question 1. Assume that po{x) > for all x E R"^. Can we say that p'^{t,x) > 
for all X G K.'' and t E [0, Tmax) ? If not, what is the maximal interval allowed for t? 

We will also recall that despite the appearance, the presence of vacuum (existence 
of zeroes of "0^) is merely a technical problem: the Madelung transform ceases to 
make sense, but a rigorous WKB analysis is available, regardless of the presence of 
vacuum. See W2.3\ 

Question 2. Let po>o e C°°(R'') with po,V0o e H''{R'') for some s > d/2 + 1. 
Suppose that the solution {v,p) to (j2.4p satisfies p{t,x) > for {t,x) G [0,t] x M^. 
Can we construct a solution to (|l.ip in (7([0,r*];/f ) (possibly with < <t)? 

We will see that in general, the answer for this question is no. Even though from 
the answer to the first question, makes sense formally, the analytical properties 
associated to (jl.ip are not as favorable as for ()2.4p . Typically, the right-hand side 
of the equation for the quantum velocity need not belong to L^(R'^). 

Question 3. Let po,(j)o G C°^(M'*) with po,V4>o G X^(M'*) for some s > d/2 + 1. 
Suppose that the solution {v,p) to (j2.4p satisfies p{t,x) > for {t,x) G [0,t] x R''. 
Can we construct a solution to (|l.ip in C{[0,t^]; X'^) (possibly with < < t)? 

2.2. Proof of Theorem 12. 3i We shall simply give the main ideas of the proof of 
Theorem 12.31 Complete proofs can be found in [35| for the cubic defocusing NLS 
equation in Sobolev spaces, and in [T] for the Gross-Pitaevskii equation in Zhidkov 
spaces. 

In the framework of this paper, we have 



Introduce formally the auxiliary function a ~ This nonlinear change of variable 
makes (j2.4p hyperbolic symmetric: 




V/(p) = Vp. 



(2.5) 




This system is of the form 




V « / 
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and the matrices Aj are symmetrized by the constant multiplier 

Id 
4 

Standard analysis (see e.g. HHj) shows that (j2.5p has a unique maximal solution 
{v,a) in C([0,rniax); Z'^), provided that s > d/2 + 1 and that, in addition, 

Tmax < +00 =^ / \\{v, a){t)\\iYl,^ dt = +00. 

Jo 

We can then define p by the linear equation 

dtp + div (pv) ^0 ; p\t=o = Po- 

By uniqueness for this linear equation, p ~ a? {a is real- valued, so p is non- negative) , 
and (u, /o) solves (j2.4p . 

We now briefly explain why compactly supported initial data lead to the forma- 
tion of singularities in finite time. The first remark is that in this case, the solution 
to (|2.4p has a finite speed of propagation, which turns out to be zero: so long as 
(w, p) is smooth, it remains supported in the same compact as its initial data. To see 
this, consider the auxiliary system p.S^ : the first equation is a Burgers' equation 
with source term 2aVa; the second equation is an ordinary differential equation 
along the trajectories of the particles. Define the trajectory by 

(2.6) ^x{t,y) =v{t,x{t,y)) ; x(0,y)=y. 

For ^ t < Tmax, this is a global diffcomorphism of W^, as shown by the equation 

^Vyx{t,y) ^Vv{t,xit,y))yyx{t,y) ; Vyx(0, y) = Id, 
and Gronwall lemma. Therefore, for a smooth function /, 

{dtf + v-\/f) {t, x{t, y)) - dt if {t, x{t, y))) , 

and (|2.5p can be viewed as a system of ordinary differential equations. 

Once the non-propagation of the support of smooth solutions is established, the 
end of the proof relies on a virial computation (like in [351 [3D], see also [H]). This 
computation shows that global in time smooth solutions to (|2.5p are dispersive 
(see also [H]). This is incompatible with the zero propagation speed of smooth 
compactly supported solutions. Therefore, singularities have to appear in finite 
time. 

2.3. A review of WKB analysis associated to ()2.ip . We consider initial data 
which are a little more general than in (|2.ip . namely 

(2.7) ^^(0,a;) = ag(a;)e''^°(^)/^ 

where the initial amplitude Cg is assumed to be smooth, complex-valued^ and possi- 
bly depending on e. Typically, we assume that there exist oq, ai € Z°° independent 
of £ such that 

(2.8) al = ao + eai+0{£^)m Z\ Ms > d/2. 
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2.3.1. First order approximation. Introduce the solution to the quasilinear system 

{dtip+hvip\^ + f{\a\^)=0 ; V'|t=o = ^o, 

(2.9) <^ 2 ^ 

1^ dta + S/ip ■ Va + —aA(p = ; a^f^Q = an- 

Theorem 12.31 shows that (|2.9p has a unique, smooth solution with a, V(/3 G Z°°. 
The main remark consists in noticing that (|2.9p implies that (V(/J, |ap) has to solve 
p.4p (oo may be complex- valued) : Theorem 12.31 viclds w,p G .2°°. We can then 
define a as the solution to the linear transport equation 

dta + V ■ Va + 2*^^ ■ v = ; a\t=o = ^o- 

Now |ap and p solve the same linear transport equation, with the same initial data, 
hence p = |ap. Using this information in the equation for the velocity, define 

fit) = ^0-1^ Qk(T)P + /(|a(r)p))dr. 

We easily check that dt (V<y9 ~ v) = Vdtf — dtv = 0, and that (cp, a) solves ()2.9p . 
Introduce the solution to the linearization of (j2.9p . with an extra source term: 

at^«+V(^-V(^(^)+2Re(aa«)/'(|ap) =0 ; f\llo = ^^ 
ata(i)+V</j-Va(i)+V(^(i)-Va+ia(i'A^+iaA(^(i) = ^Aa ; a\ll^ = ai. 

We also check that it has a unique smooth solution, with a^^\\/ip'^^^ G Z°°. The 
main result that wc will invoke is the following: 

Proposition 2.4. Let aQ,ao,ai,(pQ be smooth, wii/i ag, ao, ai, Vi^sq G . Assume 
that (EH) holds. Then 

(2.10) tA" = (ae**^*" + C (e)) e'"^/" in L°° {[0,t]; Z'), Vt < T^ax andVs ^ 0. 

This result was established in [21] when Z" ~ H''{R'^), and in [T] when Z'* = 
X'*(]R''). Note the shift between the order of the approximation between the initial 
data (known up to ©(e^)) and the approximation (of order 0(e) only): this is due 
to the fact that we consider a regime which is super-critical as far as WKB analysis 
is concerned (see e.g. [H]). In particular, the phase modulation tp^^'> is a function 
of (pojOo and ai. It is non-trivial in general, and since we are interested here in 
real- valued ao, we shall merely mention two cases (see [l6l pp. 69-70]): 

• If oi 7^ is real-valued, then ip'-^^ 7^ in general. 

• If ai = (or more generally if oi G jffi.), then ip^^^^ = 0. 

To see the first point, it suffices to notice that the equation for <^(^^ gives (recall 
that /' = 1): 

^tflHo = -2aoai. 

This shows that for (|2.2p - (|l.ip to yield a relevant description of the solution to 
(|2.ip . we have to assume ao = y/po and oi = 0. Otherwise, a phase modulation 
is necessary to describe at leading order, by (|2.10p . which is incompatible with 
the form p.2p . unless the Madelung phase (j)^ admits a corrector of order e. But 
formal asymptotics in (jl.ip give = v + 0{£'^) = V(p + 0{e'^). Hence the Madelung 
transform has a chance to give a relevant result only if ai ~ 0. 
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To check the second point of the above assertion, we set a = Re (aa'^'). Direct 
computations show that {(p^^\ a) solves, as soon as oq G M and ai G iM: 

+ • V</j(i) + 2a = ; = 0, 



9tQ; + V^s • Va = -^div (^|apV(/7(^') - aA(y9 ; an 



|t=o 

=0 = 0. 



(2.11) 



This is a hnear, homogeneous system, with zero initial data, so its solution is 
identically zero. 

To conclude this paragraph, we briefly outline the proof of Proposition 12.41 The 
approach in [Ij is the same as in j31j . with slightly different estimates. For simplicity, 
and in view of the above discussion, we assume ag = ao independent of e. We write 
the solution ip^ as = a^e^^ I' (exact formula), where we impose 

1 ^ e 

Note that both ip'^ and a*^ depend on e, because the right-hand side of the equation 
for a!^ depends on e, and because of the coupling between the two equations. Note 
also that with this approach, one abandons the possibility of considering a real- 
valued amplitude a'^. 

It is not hard to construct a solution to ()2.1ip in Z*, for s > d/2 + 2, and then 
check that asymptotic expansions are available in Z''' : 

Back to ■(/;", this yields ProDOsition l2.4l We see that the general loss in the precision 
(from 0{e'^) in the initial data to C(e) in the approximation for t > 0) is due to 
the division of by er. Note finally that even though oi = 0, one has a^^^ ^ 0: the 
corrector a^^-* solves a linear equation, with a purely imaginary (non trivial) source 
term, and so a^^-* G iM is not trivial, while ip'-'^^ = since Re (aa^^-*) = 0. 

2.3.2. Higher order approximation and formal link with quantum hydrodynamics. 
One can actually consider an asymptotic expansion to arbitrary order, 

= a + ea^i) + . . . + e^aW -f O (e^+i) , 

p' = ^ + + . . . + e"(^(") + O (e^+i) , ViV G N. 

For j ^ 1, the coefficient a^^-* is given by a linear system for {tp^^\ a'--'''), with source 
terms involving ((/j^^^ a*^'^^)o^fe^j-i . In the case oq G M and oi = 0, we know that 
ip'^^'' = 0, and a*-^^ is given by 



al'L = 0. 



9ta(2)+V^-Va(2) + ia(2)A^ + V^(2).y^^l^^^(2)^^^^(i) ^ ^^^^^ 

We check that in the case = ao G M (which includes the case where Madelung 
transform is used, ao ^ 0), all the profiles a and a^^^\ j ^ ^, are real- valued, while 
Q,(2j+i)^ j ^ g^j-g purely imaginary. Moreover, = for all j G N. This 

is formally in agreement with (|l.ip : indeed, (|l.ip suggests that cf)'^ and have 
asymptotic expansions of the form 

(2.12) + _^^2j^(2j)_^^ _ . p^~p + e2p(2)^ _^^^2j^(2,)^^^^ 
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On the other hand, we have 

(1) ^ 



« (a + . . . + £2ia(2j) + ...)'_ (ea(i) + . . . + e2j+ia(2,+i) + 

since the a(2J+i)'s are purely imaginary. This is in agreement with the second 
formal asymptotics in ()2.12p . We can check similarly that p.l2p is in agreement 
with the higher order generalization of p.lOp . in view of the special properties of 
the (f^^^^s and a^^^'s pointed out above. 

2.4. Absence of vacuum before shocks. 

Lemma 2.5. In Theorem \2.3\ assume that po{x) > for all x G (absence of 
vacuum). Then p > on [0,T,„ax) x K''- 

Proof. As in Section [2?2l we use the fact that on [0,Ti„ax), the equation for the 
density is just an ordinary differential equation. Introduce the Jacobi determinant 

My) = dcWyx{t,y), 

where x{t, y) is given by (|2.6p . We have seen that Jt{y) > for {t, y) G [0, Tmax) x 
R'^. Change the unknown p to r, with 

r{t,y) = p{t,x{t,y)) Jt[y). 

Then for ^ i < Tmax, the continuity equation is equivalent to: dtr = 0. Therefore, 

^) = T / ], nn Po {y{t. x)) , 
Jt {y{t,x)) 

where x i— ?> y(t, x) denotes the inverse mapping of y i—> y). □ 
Wc infer: 

Proposition 2.6. Under the assumptions of Theorem \2.3[ assume that po{x) > 
for all X G (absence of vacuum). Let < T < Ty^ax, o-nd K be a compact set in 
R'^. There exists e{T, K) > Q such that forQ<£^ e{T, K), > on [0, T] x K. 

Proof. Proposition 12.41 shows that 

IV/ 1 = \a\+0{e) in L°° ([0,T] x M^) . 

Note that the constant involved in this 0{e) depends on T in general. Recalling 
that a = y^, Lemma 12.51 shows that 

min a(i, x) — c(T, K) > 0. 

it,x)e[0,T]xK 

Now for < £ ^ £(r, K) < 1, 

M ^ lc{T,K) in ([0,T] x R'') , 

and the result follows. □ 

In the case a; G T'', this shows that before the formation of shocks in the Euler 
equation, and provided that e is sufficiently small, the amplitude remains positive: 
the right-hand side of (|l.ip makes sense. This point was remarked initially in 
P7] . Note that the result of [55] in the one-dimensional case x G [0, 1] shows that 
suitable boundary conditions lead to the existence of finite time blow-up for (jl.ip . 
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Therefore, the above result is quahtatively sharp (quahtatively only, for it might 
happen that the solution to (jl.ip remains smooth longer than the solution to (|2.4p ). 

Finally, we show that the compactness assumption in Proposition 12.61 can be 
removed in the case of the Gross-Pitaevskii equation. As regards the nonlinear 
Schrodinger equation on K'', this issue seems much more delicate and will not be 
addressed in this paper. Assume that the Gross-Pitaevskii equation is associated 
with the boundary condition at infinity 

\^Ht,x)\ 1. 

|a;|— foo 

Such a condition is used frequently in physics, possibly with a stronger one, of the 
form (see e.g. [35] and references therein) 

>-oo 

for some fixed asymptotic "velocity" v°° £ Mf^. 

Now, as regards the Gross-Pitaevskii equation, putting together the continuity of 
tp'^ over [0,T] x M."^, the compactness of the time interval [0,r] and Proposition [mi 
we get 

Corollary 2.7. Under the assumptions of Theorem \ 2.3l assume that pa{x) > 
for all X e (absence of vacuum). Assume moreover that the Gross-Pitaevskii 
equation is associated with the boundary condition at infinit^ 

\rit,x)\ 1. 

\x\^OQ 

Let <T < Tniax- There exist e(T),c(T) > such that 

^ c(T) on [0, T] X W^, for all 0<e ^ e{T). 

2.5. Functional spaces associated to the Madelung transform. It is rather 
easy to see that the answer to Question [5] is no, in general. Consider for po the 
function in the Schwartz class 

Po(x) = e-l^l", fc>l. 

At time t ~ 0, the quantum pressure (right-hand side of grows like |a;p''~^ 

hence the velocity v'^ has no chance to belong to H" for general initial data in H°° . 
Thus, working in Sobolev spaces for general initial data does not make sense for 
()l.ip , while the results in [31] show that it is a fairly reasonable framework to study 
the semi-classical limit of p.ip . 

On the other hand, like for the absence of vacuum, the answer to Question [3] is 
positive, at least if we consider some special boundary conditions at infinity. 

3. Solving the QHD system by a direct approach 

In the present section, we describe an efficient method to solve directly the 
hydrodynamic form of (|2.ip given by (|l.ip , once performed the Madelung transform. 
This method enables us to study the corresponding initial value problem for (|l.ip 
with data {vo,po) such that (?;o,V/5o) has a high order Sobolev regularity and pQ 
is positive and bounded away from zero. In addition to local-in-time well-posedness 
results, we get (see Theorems 13.11 and 13.31 below) nontrivial lower bounds on the 



Note that this imphcs that there exists c > such that poi^) ^ c for all x £ ' 
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first appearance of a zero for tlie solution, wliicli are of particular interest for the 
study of long-wavelength asymptotics if the data are a perturbation of a constant 
state of modulus one. 

Wc here closely follow the approach that has been initiated in [8]. To help the 
reader to compare the present results with those of the previous section however, 
we keep on using the semi-classical scaling given by ()2.1|) (whereas e = 1 in [8]). 

The use of a suitable extended formulation for and of weighted Sobolcv 

estimates will be the key to our approach. Let us stress that, recently, similar ex- 
tended formulations have proved to be efBcient in other contexts for both numerical 
(see [21]) and theoretical purposes. As a matter of fact, in the last paragraph of 
this section, we shall briefly explain how our approach based on such an extended 
formulation carries over to the more complicated case of Korteweg fluids. 

3.1. Solving the QHD system by means of an extended formulation. The 

"improved" WKB method that has been described in the previous section amounts 
to writing the sought solution ip'^ as: 

for some complex valued function and real valued function ip"^ . In this section, 
we rather start from the Madelung transform 

i/j" = y/p^e^'^^ where Vp^ = IV""!, 

then write 

= e'*^/' with = 0^ - i| logp*^, 

and consider the redundant system that is satisfied by both and — V^"^ — 
-\- iw^ with i;^ = V(/)^ and ~ —■^Vp'^. 

In order to obtain the system for z"^ , we first difFcrcntiatc the density equation 
in (ITTTI) . This yields 

dtw^ -I- V(w^ • w") = Ivdivw". 



Next, we notice that 



■ div ^ — \w^ 



2 y/y 2 2' 

In consequence, the equation for tj^ rewrites 

dtv" + v" ■ S/v" ~ ^\/\w'f + |V div w' + Vfip") = 0. 

Of course, z^ is a potential vector-field, hence V div z^ ~ Az'^ so that we eventually 
get the following "extended" system for (p'^,z'^): 

(3 1) I dtzJ + iv(z^' • zJ) + Wfip^ = ^|Az^ 

1 dtp"" + dW{p''v^) = 0, 

d 

where we agree that a ■ b ajbj for a and 6 in C*. 
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Let US now explain how Sobolcv estimates may be derived for (p^, z^) in the case 
where = 1 + 6^ for some 5^ going to at infinity, and /(p^) = (to simphfy). 
The following computations are borrowed from [S]. For notational simplicity, we 
omit the superscripts e. 

In order to get the basic energy estimate, we compute^: 

/ ((l + 6)|z|2 + 52) =2f / {l + b){z,dtz)+ ( bdtb)+[ dtb\z\^ 

h h h 

d 

where the notation (a, b) = Re aj Re bj + Im Oj Im bj has been used in Ii . 

Further computations yield Ii = Ii i + 112 + ^1,3 + /i,4 + ^1,5 and I2 ~ I2.1 + ^2,2 
with 

11,1 = - l{z,Wb), h.i^- I bdivv, 



2 



b{z,\/b), I2.2 j ^ div(6z;). 



A,3 = y (2:,i-Az), 
A,4 = y 5(z, j|Az), 
A,5 = / p(z,V(z-z)), 



2 _ 

Using obvious integrations by parts we readily get 

-^1,1 + -^2,1 = 0, /i,2 + ^2,2 = and /i,3 = 0. 
Therefore, integrating by parts in /i^4 also, we get 

^ y"((l + + b^)=- J ^{z, *V2 ■Vb)-2 j p{z, Vz-z)+ J dtb 

For "general" functions b and z, the appearance of the terms V6 and Vz would 
preclude any attempt to "close" the estimates. In our case however, as the algebraic 
relation — eV6 = 2pw holds true, one may avoid this loss of one derivative for one 
may write 

-e{z,iVz ■ V6) - 2p(z, Vz • z) = -2p(z, Vz • v). 
Now, integrating by parts an ultimate time, we conclude that 

I J{p\z\' + b')=-JpvV\z\' + |a,5|zp = o. 

Hence, J (p|zp + 6^) is a conserved quantity. 

The same algebraic cancellations may be used for getting higher order Sobolev 
(or Besov) estimates. Indeed consider an "abstract" pseudo-differential operator 



The method may seem uselessly complicated. However, the algebraic cancellations that are 
going to be used remain the same when estimating higher order Sobolcv norms. 
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A{D) (for instance a differential, a fractional derivatives or a spectral localization 
operator). Then one may write 

I l{{l + b)\AiD)z\^ + {AiD)bf) 

= 2(j{\ + h){A{D)z,dtA{D)z)+ 1 A{D)bdtA{D)b\ + J dtb\A{D)z\\ 



ii I2 I3 

We notice that Ii = + Ii^2 + h.z + Iia + and I2 = h,! + ha with 

Ii.i= - j {A{D)z,VA{D)b), /2a=-y" A{D)hdiYA{D)v, 

h.2 = ' J h{A{D)z, VA{D)b), l2,2 = - J A{D)bAiD) div(&t;). 
/i.3= J {A{D)z,t^AA{D)z), 
/i,4= J b{A{D)z/-AA{D)z), 

h.B = -\ I p{A{D)z,A{D)W{z-z)), 

As above, obvious integrations by parts ensure that /i.i + l2,i = and /i_3 = 0. 
Next, using again integrations by parts, we notice that 

h,2 + h,2 = j A{D)b div[6, A{D)]v. 

Finally, integrating by parts in /1.4 and using the fact that — ffV6 = 2pw yields 

h.A = j {A{D)z,i{yA{D)z) ■ [pw)), 

and we have 

/i,5 = - J{A{D)z, A{D)\/ z ■ {pv + ipw)) + Jp{A{D)z,\/A{D)z ■ z- A{D){\/z ■ z)). 
Therefore, using the fact that dtp + div(pf ) = 0, we conclude that 

2(/i,4 + /i,5) + ^3 = 2 y" p{AiD)z, yA{D)z ■ z - A{D){Vz ■ z)). 
Putting all the above equalities together, we thus get 
~ J{p\AiD)z\' + iAiDW)= I A{D)bdW[b,AiD)]v 

+ I p{A[D)z,VA{D)z-z-A{D){\/z-z)). 



If, say, A{D) is a fractional derivatives operator, then one may show by means of 
classical commutator estimates that the right-hand side may be bounded by 

CWpWl^ II [Db, Dz)\\l^ \\{A{D)b, A{D)z)\\l. . 

Therefore, 
d 



dt 



{p\A{D)z\^ + {A{D)bf) ^ C\\p\\L^\\{Db,Dz)U^\\{A{D)h,A{D)z)\\ 



2 

L2- 
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Denoting by (t) the left-hand side and resorting to GronwaU lemma, we thus get 

(3.2) EAit) ^Ea{0)c^p(c \\p\\L^\\p-'\\L^^\\{Db,Dz)\\L^.dTy 

It is now clear that whenever Db and Dz arc bounded in L^{[0,T]; L°°) and p is 
bounded from below and from above then we get a control of A{D)b and A{D)z in 
L-([0,r];L2). 

Taking A{D) = A" and performing a time integration, (|3.2[) implies that 

(3.3) \\{b,z){t)\\H~.<c(\\{bo,zo)\\H^+ r \\{Db,Dz)\\L^\\{b,z){t)\\Hsdt' 



(3.4) 



for some constant C = C{s,d, \\p^^\\l°^). 

So assuming that s > 1 + d/2 and using Sobolev embedding and Gronwall's 
inequality, one may conclude by elementary methods to the following statement. 

Theorem 3.1. Let s > I + d/2. Assume that po = l + ba for some bo e H''+'^{R'^) 
such that 1 + &o > 0, and that vq £ H''{M.'^). Then there exists a time 

T^To:= II With C^C{s,d,\\p^^L^) 

such that p.ip has a unique solution {v^ , p^) on [0, T] x with p~ = 1 + b^ bounded 
away from and (v^6'') € C([0,r];i?" x iJ^+i) n C^i[0,T]; H''-^ x i7"-i). 

Remark 3.2. Combining basic energy estimates for the wave equation with the 
above result, one may control the discrepancy between (6, v) and the solution to 
the acoustic wave equation 

( dtv + Wb = ; vit=o = vo, 

[ dtb + divv^Q ; fe|t=o = ^o- 
We have, up to time Tq, 

|j(i;^'-i.,6^-6)(t)||H.-2 <C(i|l(5o,^xo)!pH.+ixif-' +e'^ll(^o,^io)llH=+ixHO- 
Note also that (|3.3p provides a blow-up criterion involving the W'^'°° norm of 
b and the Lipschitz norm of v. In particular, this implies that for given data in 
H'^ (s > 1 + d/2), the lifespan in i/* is the same as the lifespan in W , for any 
I + d/2 < s' < s. 

3.2. Dispersive properties and improved lower bounds for the lifespan. 

The system for (6, v) reads 

(3.5) ( 5*- + -V. + V6=^v(iA6-^|V6P 
I dfb + div V = — div(5w). 

Therefore the linearized system about (0, 0) is not (|3.4p but rather 

(3.6) I dd' + m^'-^S/Ab, 

[ dtb + div 0. 

A straightforward spectral analysis (based on the Fourier transform) shows that 
the above linear system behaves as the wave equation with speed 1 for frequencies 
small with respect to 1/e, and as the Schrodinger equation with coefficient e/2 in 
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the high frequency regime. In fact, in dimension d 2, it is possible to prove 
Strichartz inequahties (related to the wave and Schrodinger equations for low and 
high frequencies, respectively) for (|3.6I) . In the case of small data {bo,vo), these 
inequalities allow to improve the lower bound for the lifespan (see also |22j where a 
similar idea has been used in the context of the incompressible limit for compressible 
flows). For the sake of simplicity, let us just state the result in dimension d ^ 4 
(the reader is referred to for the case d = 2,3 and for more details concerning 
the approximation of the solution by p.6p ): 

Theorem 3.3. Under the assumptions of Theorem \3.1\ with s > 2 + d/2, then the 
lifespan T may he hounded from helow by 



may be bounded in terms of t and of the data, up to time Ti . 

3.3. Extended formulation for Korteweg fluids. Compared to the "improved" 
WKB method, the main drawback of the direct approach based on an extended 
formulation for solving p.ip is that vanishing solutions cannot be handled. 

On the other hand, the direct method is robust enough so as to be used to 
solve locally more complicated models such as the following system governing the 
evolution of inviscid capillary fluids: 



[ dtp + div{pv) = 0. 

Physically, the function n correspond to the capillary coefficient. Obviously, Sys- 
tem (jl.ip is included in p.7p (take k{p) = £^/(4p)). In the general case, introducing 



Note also that in the general case, the second order term \7{a{p) div z) is degenerate. 

The case of System (jl.ip is particularly simple inasmuch as a is the constant 
function e/2 and V divz = Az. 

For general capillarity coefficients, one may prove a local well-posedness result, 
similar to that of Theorem 13.11 This has been done in [7] . The proof relies on the 
use of weighted Sobolev estimates, with a weight depending both on p and on the 
order of differentiation. 

As for the QHD system, we expect the potential part of the solution to Sys- 
tem (j3.7p to have dispersive properties in dimension d ^ 2. The general situation 
is much more complicated however, because those properties are related to those 
of the quasilinear Schrodinger equation. To our knowledge, this aspect has been 
investigated only very recently in a work by C. Audiard [5] that concerns potential 
flows. 




(3.7) 






iV(z.z). 
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4. ASYMPTOTICS FOR THE GrOSS-PiTAEVSKII EQUATION 

This section is concerned with the existence and asymptotics of travehng wave 
solutions for the Gross-Pitaevskii equation 

(4.1) zVt + AV^ + (1 - iV^nv- - 0, 

which may be obtained from (|2.1|) (with /(r) = r — 1), up to the factor after 
performing the change of unknown: 

Equation (|4.ip is associated to the Ginzburg- Landau energy (or Hamiltonian) : 

(4.2) nw^l I ivv'p + i / (i-l^n^ 

As a consequence, in contrast with the cubic NLS equation, the natural energy 
space for is not ^^[W^) but rather 

E{M.'^) = {i/- e i?ioc(K''), s.t. •H(V') < +oo}- 

As pointed out before, for H{'ip) to be finite, |-!/'| must, in some sense, tend to 
1 at infinity. This "nontrivial" boundary condition provides (j4.ip with a richer 
dynamics than in the case of null condition at infinity which, for a defocusing NLS 
type equation, is essentially governed by dispersion and scattering. For instance, 
in nonlinear optics, the "dark solitons" are localized nonlinear waves (or "holes" ) 
which exist on a stable continuous wave background. The boundary condition 
|V'(i,a;)| — >■ 1 at infinity is due to this nonzero background. In the context of 
superfluids, 1 is the density of the fluid at infinity. 

Similarly to the energy, the momentum 

is formally conserved. This quantity is well defined for -0 g H^iW^) but not for 
solutions with a finite Ginzburg-Landau energy. A major difficulty in the theory of 
the Gross-Pitaevskii equation is to find an appropriate definition of the momentum 
which leads to a conserved quantity. A natural definition would be 

but this would require for instance that ?/' — 1 e L'^{W^). 

In any case we will denote by p the (scalar) first component of V which will play 
an important role in this section. 

Recall that the use of the Madelung transfornfl 

i> = ee'\ 

leads to the following hydrodynamic form of the equation for q and v — 2V0 : 

{atw + w • Vw + 2Vq^ = 2V ( — I , 
V ^? / 



'in this section, we use the normahzation of [9lll0|, instead of tp = y^e"^. 
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As pointed out before, if we discard the right-hand side of the first equation and 
fook at as the density of a fluid with vefocity w, then the above system coincides 
with the Euler equations for a compressible fluid with pressure law P{p) = p^. In 
particular, the speed of sound waves near the constant solution i; = 1 is given by 

As wc will see below, this sound speed (the value and relevance of which is not so 
obvious if looking at the original Gross-Pitaevskii equation (|4.ip ). plays an impor- 
tant role in various aspects of the dynamics of (j4.ip . The value of Cg may be also 
found by neglecting the quantum pressure term and linearizing for a perturbation 
^ = (1 + ^) exp(i(/)). This leads to the wave equation: 

Note that if the quantum pressure is included (as in ()3.6p ) then the linearization 
reads: 

dfg-2Ag^A'^g = 0, 
which, roughly, is the factorization of two linear Schrodinger operators. 

The Madelung transform and the hydrodynamic form of (|4.ip turn out to be of 
great interest to study the Gross-Pitaevskii equation with finite Ginzburg-Landau 
energy since the solution is expected to have very few "vortices" (or cancellations). 
Even in the "Euler limit" that has been presented in Section [21 one can use it 
outside the vortices ([Mj [15]) to study the traveling waves of sufficiently small 
velocities. Let us also stress that the hydrodynamic form of the (one dimensional) 
Gross-Pitaevskii equation is needed in order to define a generalized momentum in 
the context of the orbital stability of the black solitons (such solitary waves have 
zeroes. . . ); see [IT] . 

In the present section, we shall concentrate on the transonic limit of solutions 
to the Gross-Pitaevskii equation. We shall first present a result pertaining to the 
asymptotics of traveling waves with speed c tending to the sound speed , in the 
case d = 2, in connexion with the (KP I) equation (see below). Next, for the one- 
dimensional case, we give an accurate description of the transonic long wave limit 
of (|4.ip in terms of solutions to the KdV equation. 

4.1. The transonic limit of finite energy traveling waves. Finite energy trav- 
eling wave solutions of (j4.ip are solutions of the form ip{x, t) ~ 9{xi — ct, a;-*-) where 
7i{9) < +00 and denotes the transverse variables X2,--- ,Xd- The profile 9 
satisfies the following equation: 

(4.4) icdi9 + A9 + 9{1 - \9f) = 0. 

A suitable functional setting for the study of such traveling waves is the space: 

W{M.'^) = {1} + V{R'^), 

with 

V{R'^) ={i^:R''^C, (VV',Re?^)€L2(]R^)^Im?/'ei''(M''), VReV' G L^/^{m.'^)}- 

Indeed, given that W{R'^) is a subset of the energy space E{M.'^), for any data in 
W{W'-), Equation (|4.ip admits a unique solution. In addition, one may show that 
this solution stays in M^(R'^) (see [28l[29]). Furthermore, the quantity {iditp, ■0—1) is 
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iiitegrable whenever ip S W{W^), so that the scalar momentum p(V') is well-defined. 
This is a consequence of the identity 

(4.5) {idiil),il)-l) =ai(RcV')ImV'-5i(Im7/')(RcV- 1), 

and various Holder's inequalities. 

So finally, T-L and p are continuous on W^(IR'^) and all finite energy subsonic 
solutions to (|44)) have to belong to WiW'-). Moreover, if V' G VF(M'') may be lifted 
as -0 = £iexp(i(/)) then 

(4.6) pW = \f > - 1) = i / (i-e')ai0. 

Notice that for maps which may be lifted, with g ^ ^, the last integral makes sense, 
even if we assume that tp only belongs to the energy space E{R'^). 

To simplify the presentation, we shall focus on the simpler case d = 2. We shall 
denote x = xi and y = x^ = X2- 

It is proven in [S] that, for any p > 0, the following minimization problem 

(4.7) HnuniP) = inf{H(VOW^ e W{R^),p{i;) = p}, 

has a solution Up which is a nontrivial traveling wave. We call it a ground state. 

In the rest of this subsection, we focus on the asymptotics p going to for Up, 
in connexion with the Kadomtsev-Petviashvili I (KP I) equation 

ut + uux + Uxxx - d^^Uyy = 0, 

where the antiderivative is defined in Fourier variables by dx^f{£,) = j^fiO- 

The following proposition states that the corresponding speed c(up ) tends to Cs 
and gives the first term in the asymptotic expansion for both 'Hmin(p) and c{up). 

Proposition 4.1. There exist positive constants pi, A'g and Ki such that we have 
the asymptotic behaviors 

(4.8) ^^p-^ - AV sS \/2p - H,„i„(p) ^ AV, VO^p^pi, 

^KP 

where Skp stands for the action of the KP I ground state N of velocity 1, that is 

(4.9) Skp^\I {dxNf + \l {d-\dyN)f^\f N^ + \! N\ 

Moreover, the map Up has no zeroes and there exist some positive constants p2, 
K2 and A'3 such that 

(4.10) A2p2 y2-c(wp) < A3P^ V0s;p<p2. 

Actually, it was established in [5] that if ij.' is a finite energy traveling wave of 
sufficiently small energy, then 

^<I^I<1- 

This implies that small energy traveling waves have no zeroes and thus can be 
lifted according to the Madelung transformation. This is in particular the case of 
minimizers Up corresponding to small enough values of p. 
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The transonic limit to KP I for travehng waves is obtained through the fohowing 
change of scales: 



We then set 



i = eWa;, y = -^y- 



1 I / i2 AT f \ 6 I ^ 



It turns out that N^i, converges to a traveling wave solution of the (differentiated) 
KP I equation as £(-0) —J- 0, that is it approximately solves the equation 

(4.11) ~Wxx-Wyy+Wxxxx + {wWx)x^Q- 

More precisely, denoting Np = 1 — |wpp where Up = |up|e*'^p is a minimizer of the 
energy with fixed momentum p, and 

n ( \ ^^A. ( ^ 

we have (see [lO]): 

Theorem 4.2. There exist a subsequence (p,i)n6N tending to 0, as n ^ +oo, and 
a ground state w of the KP I equation such that both Np^ and Qp^ tend to w in 
p^fc,gj-]|j2-j ^j^^ If ^fq fj^^fi q g (X, +oo]) as n goes to +oo. 

Remark 4.3. Those results on the transonic limit of solitary waves have been 
recently extended in [19] to the three-dimensional case with also general nonlinear- 
ities (see also [S^ for other existence results). The additional serious difficulty is 
that the Gross-Pitaevskii ground states solutions are no longer global minimizers 
of the energy with fixed momentum and thus not expected to be stable. 

4.2. The unsteady transonic long wave limit of the Gross— Pitaevskii equa- 
tion. The Madelung transform is also crucial to derive and justify the transonic 
(weak amplitude, long wave) limit of the Gross-Pitaevskii equation. 

With the scaling which is used in this section and for data which are perturba- 
tions of order e'^ of a constant state with modulus 1, Theorem 13 . 1 1 and the remark 
that follows ensure that the linear wave equation gives a good approximation of the 
solution for t = o{e^^). In the present paragraph, we describe what happens at next 
order. We shall see in particular that, up to times of order e^^, the Korteweg-de 
Vries (KdV) equation 

Ut + UUx + Uxxx = 0, 

gives an accurate approximation of the one-dimensional Gross-Pitacvskii equation 

(4.12) idt^ + dxx^ + (l ~ = 

with data which are small long-wave perturbations of the constant one, namely 
ijj = ge"<' with 

/ 2 \ 1/2 

(4.13) go=i^l- jN^iex)j , 00 = ^©"M- 



^Actually, to a ground state solution of the KP I equation, that is a minimizer of the Hamil- 
tonian with fixed norm. 
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for < £ <C 1, and and = OxQ^ uniformly bounded in some Sobolev 
space iJ'^(M) for sufficiently large k. 

Wc rcfcr to [TH [13] for a detailed analysis and will only summarize the limit to 
(long) waves propagating in two directions and following a coupled system of KdV 
equations. 

Recall that the one-dimensional Gross-Pitaevskii equation ()4.12|) is globally well 
posed in the Zhidkov type spaces Y'^, 

Y''{R) = {i,eLl^{R;C), 1 - e ^^(M), d^i; e H''-\R)} , 

for any integer k ^ 1. 

Moreover the Ginzburg-Landau energy 'H{ipit)) is conserved by the flow and, 
provided 'H('0o) < the corresponding solution does not vanish so that one 
may write ip — giexp(i0), for some continuous function (f). 

We consider data as in (|4.13p with small enough e and assume in addition that 



m\Mi 



\dxe°\\Mm < +0O. 



Here, || • \\m{R) denotes the norm defined on Ll^^{R) by 



l/ll 



M( 



sup 

{a,b)eM 



f{x)dx 



We next introduce the slow coordinates 

^ e{x + V2t), a;+ = e{x - \/2t), and 



2V2 



The definition of the coordinates x~ and x^ corresponds to reference frames trav- 
eling to the left and to the right, respectively, with speed \/2 in the original coor- 
dinates {t,x). We define accordingly the rescaled functions and <df as follows: 



(4.14) 

Setting 
(4.15) 



N^{t,x^) 



6 



r]{t,x) 



6 / 2^2' 



T X^At 



with rj = 1 — , 



ef{r,x^) = ^<j>{t,x) 



6^/2 / 2V2T ^ _^ 4t 



C/-(t,.t-) = - {N-iT,x-)+dx-e-iT,x-)) , 
U+{t,x+) = i {N+{t,x+) - a,^e+(r,x+)) , 



the main result is (see [T3] for details): 

Theorem 4.4. Let k ^ and e > be given. Assume that the initial data ipo 
belongs to y'^+^(R) and satisfies the assumption 

\\N°\\Mim + \\d.Q"AMim + \\N^\\H^+^m +e\\d',+'^N^\\mR) + Wd^QtWH-^nm) < Ko. 

Let and denote the solutions to the Korteweg-de Vries equations 
drU- + dlM- +U-dx-U- = 0, 

and 

drU+ - dl+U+ - U+dx+U+ = 0, 
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with the same initial value as and respectively. Then, there exist positive 
constants ei and Ki, depending only on k and Kq, such that 

\\U-{t,-) -U-i;T)\\H^^K^ + \\U+{t,-) -U+{;t)\\h>^^W) ^ Ki£^ cxpiKi\T\) , 

for any r G M provided s ^ ei. 

We now turn to the two (or higher) dimensional case, which is studied in |20| 
for a general nonlinear Schrodinger equation of the form similar to (|4.ip 

where / is smooth and satisfies /(I) = 0, /'(I) > 0. 

One also uses a "weakly transverse transonic" scaling, namely 

T ce^t, Xi e{xi - ct), Xj = e^xj, j = 2, • • • ,d. 

After performing the ansatz 

ij^t, X) = {1+ e^A^t, X)) exp(ie(/)^(t, X)), 

the hydrodynamic reformulation of the Gross-Pitaevskii equation is used to recast 
the problem as a singular limit for an hyperbolic system in the spirit of |31] . Then 
smooth iJ* solutions are proven to exist on an interval independent of the small 
parameter e. Passing to the limit by a compactness argument yields the convergence 
of the solutions to that of the KP-I equation. Note however that this method 
does not provide a convergence rate with respect to e, contrary to the KdV case 
considered above. 

In comparison, for such data. Theorem 13.31 would ensure that the linear system 
(|3.6|) gives a good description of the solution only for times that are o(e~'^) (see the 
introduction of [5] for more details). 

5. Global existence of weak solutions to a quantum fluids system 

We aim at providing an elementary proof of the result by P. Antonelli and 
P. Marcati in [Sj [4] concerning global finite energy weak solutions to the QHD 
system. Here is the statement: 

Theorem 5.1. Let the initial data (po,Ao) G W^'-^ x be "well-prepared" in the 
sense that there exists some wave function ipQ G such that 

Po = iV'oP and Jq ^/poAo = Im(7/!oV'0o). 

Assume that f{r) ^ r"^ for some integer a such that W'^-^{W^) ^ L'^+i(R'*). 

There exist some vector-field A g L°°(R;L^) and some nonnegative function 
p e L°°{R; n 1"+^) with WJp e L°°(M; L^) such that the following system holds 
true in the distributional sensqj; 



(5.1) 



dtp + div J = 0, 

dtJ + div (A A) + V(P(p)) = ^ AVp - div (V^p ® V^p), 
djj'' - dkJ^ = 'i'ts^dj^fP - '^J^'dk^/p for all {j, fc) e {1, • • • , d}^ 

(p, J)\t=Q = (PO, ^o), 



^In the smooth non-vanishing case, the right-hand side of the second equation coincides with 
that of the velocity equation multiplied by p in HI. It with £ = 1, and the third equation just 
means that there exists some function <f) such that J = p^cft. 
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With J := and P{p) ;= pf{p) - F{p) with F{p) = f{p') dp'. 
In addition, the energy 

is conserved for all time. 

Proof. Let us first prove the statement in tlie smootli case, namely we assume tiiat 
tlie data ipQ is in for some large enough s. It is well known that 

(5.2) zatV' + ^A^. = / (|^|2) ^ ; ^|t=o = ^o 

has a unique solution ij) in C(R; H^) whenever s ^ 1 (see e.g. [TBI ISHl El] ) a-nd that 

(5.3) VteR, / (\\Vm?+F{\m?)]^ I fJlV^ol' + FdV'ol') 
Let us set p := \iJj\^ and 

|V'(a;)|-V(a;) if + 0, 

if i}{x) = 0. 



(5.4) ^(a;) 



We claim that 

(5.5) VVp = Re(<^V?/') a. e. 



Indeed, for e > 0, let us set :— i^l + Then ((/)e) converges pointwise to 
(j) and an easy computation shows that 

(5.6) (pe"^ — 5" -v/p uniformly. 
Next, we compute 

(5.7) V(^eV') = Re(0eV^/') + Rc(V'V0e)- 

The first term in the right-hand side converges pointwise to Re(0V'0) hence in L\^^ 
owing to Lcbesguc's theorem as it is bounded by IVV'I- 
As for the last term, we have 

^^i^^e) = (^,^[',p^3/2 Rc(^v^). 

It is clear that the right-hand side converges pointwise to and is bounded by \ 
Hence it also converges to in L\^^. 

So finally, putting these two results together with (|5.6p and ()5.7p . one may con- 
clude to ([53]) . 

Let A := Im(^V'(/') and J := ^/pA. We claim that (p. A, J) satisfies (|5.ip . Indeed, 
from (|5.2p . we see that 

dt{\il)\^) = 2Re(MV') = -MVJAV') = -Imdiv(V;V?/'), 

hence 

(5.8) dtp div J = 0. 
Next, we compute 

a, j'= - 9fcj^ = im(aj(v;afcV) - 5fc(MV)) = 2im(a,^afcV). 
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Recall (see e.g. Theorem 6.19 in [Sj) that 

(5.9) V?/- = a. e. on ^-"^({0}) 
whenever ip is locally in W^'^. 

Therefore, given that = 1 on -0^^(C \ {0}), one may write a. e. 

So we eventually get 

(5.10) djj'' - dk.P = 2a,VP - 2dk^pkK 
Next, we compute 

dtJ^ Im(9tV; VtA + ij\/dtip) . 
Hence, using the equation satisfied by ip, we get 

dtJ = I Rc(^VA^ - A^Vt^) + Re(/(|V'n^VV' - ^pVifi\if)iP)) . 

A B 

If ■0 is C'^ and the function / has a derivative at every point of IR+ , then straight- 
forward computations show that 

(5.11) S = -p/'(p)Vp = -V(P(p)). 
Next, we see that (still in the smooth case) 

VAI^Z-I^ = 4Re(vV : V^/J) +2Re(VV'AV') + 2Re(V5VAV'). 

Hence 

A = ivAlV-P - 2Re(A?/;VV') - 2Re(vV : V^). 

So we get 

(5.12) = ivA|7/>|2 -div Rc(VVi(g) V?/-). 

Now, using again (|5.9p . one may write at almost every point of R'', 
Re(V'0 (g) V'0) = Re((/)V'0 «) ^V?/-), 

= Re{(t)Vip) «) Re((^VV') + Im(0V'0) Im(^V?/'). 

Therefore, we have 

(5.13) Rc(V'0® V-iA) = V^p® V^ + A(g) A. 
Putting this together with (|5.1ip and (|5.12p . one may conclude that 

(5.14) dt J + div (A ® A) + V(P(p)) = ^ AVp - div (V^^^ ® V^^S). 
Of course, as owing to (j5.9p 

(5.15) IVtAI^ = (Rc(0VV'))^ + (Im(0VV))2 = iV^j^l^ + |A|2 a. e., 
the energy equality for ijj recasts in 

(5.16) £(i|A(t)|2+i|Vv/;^p+F(p(t))) =£(i|Aop+i|Vy;^P+nPo)). 
This completes the proof in the smooth case. 
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Let us now treat the rough case where i/'o belongs only to H^. Then we fix some 
sequence (V'o,n)neN of functions in H" (with s large) converging to i/'o in H^. Let 
us denote by V'n the solution of (|5.2[) in C{M.;H'^) corresponding to the data ■00, m 
and by -0 G C(M; H^) the solution to ((5?2|) with data ipo- 

From the first part of the proof, we know that there exists some sequence (0n)n6N 
of functions with modulus at most 1 such that if we set A„ :— Im(0„V0„), p„ :~ 
IV'nP and J„ := y/p^A,-, then 

VVp^ = ^^(^nVi/'n) in L^, 

and (p„, A„, J„) satisfies (|5.ip with data (po,n, >/o,n), together with the energy equal- 
ity 

(5.17) £^ (^|A„(OP + ^|Vv/;^|' + F (p„(t))) 

= (^1 Ao.nl' + i IVVp^P + F{po,„)) . 

We now have to prove the convergence of (p„, A„, J„) to some solution (p, A, J) of 
(|5.ip satisfying the energy equality. 

On the one hand, standard stability estimates (based on Strichartz inequalities) 
guarantee that 

(5.18) ^n^^ in L^,iR;H^). 

On the other hand, because (0„)„gN is bounded by 1, it converges (up to extraction) 
in L°° weak * to some function (j) such that ||(?I>|jL'° ^ 1. As V^/'n converges strongly 
to V-^ in L^, this implies that 

(5.19) 0„VV'„ in 

In turn, as obviously — > in and as V^^p^ = Re(0„V-!/'„), we deduce 
that 

V^ = Re(0VV') and V^P^ ^ V^P in 

Given that A„ = Im(0„VV'n), (|5.19p also ensures that 

A„ ^ A := Im(0V0) in L^. 

In order to establish that strong convergence in holds true, it sufiiiccs to show 
that (up to an omitted extraction) 

(5.20) WV^Wl- ^ l|Vyp||i2 and ^ I|A||l^. 
On the one hand, the weak convergence ensures that 

IIVVpIU^ 5% liminf ||V,/p;r||L2 and j| Aj|i2 ^ liminf || A„j|i2; 

on the other hand, given that (V'i/'n)rieN converges strongly to 'Vip in and that 
(|5.15p holds true for ip and ipn, we may write 

!|vVp|li^ + !|A!li2 = ||vv.||i2, 

= lim |iV0„|li2, 

n— )-+oo 

= lirn (WS/VP^Wh + WAJl,), 

> (liminf \\\7^\\L2f + (liminf ||A„||i2)^ 
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Therefore ()5.20|) is satisfied. As a conclusion, we thus have established that 

\f(hi ^ \fp in and A„ A in . 

Of course, this implies that J '■= -^/p^ iii so it is easy to pass to the limit 

in (|5.ip and in the energy equality (|5.16p . The details are left to the reader. □ 

Appendix A. Conservation laws 

In this Appendix we review conservations laws for the nonlinear Schrodinger, 
QHD and compressible Euler equations. Even though most of the results are clas- 
sical (as concerns the Schrodinger equation, they may be found in the textbooks 
[431 144] for instance; links between Schrodinger and Euler conservation laws may be 
found in |16j). we believe the relationships between the aforementioned equations 
to be of interest. In addition, those conservation laws are still meaningful for the 
less classical framework of general Korteweg fluids. 

A.l. The case of Schrodinger, QHD and compressible Euler equations. 

For the time being, we consider the following system 

[ dtp + div {pv) = ; p\t=o = po, 

which is the QHD system if e > 0, and the compressible Euler equation if e = 0, 
and the nonlinear Schrodinger equation: 

(A.2) je9t^4.^A^^/(|^|2)^ . ^-it^o = V'o- 

Recall that for e > 0, one may pass formally from (jA.2[) to (jA.ip by setting 

iP = T^e*"^/^ and v = V0. 

In what follows, the function / is assumed to be continuous on and, say, 
on (0, +oo[, standard cases being /(?-) = t"' and /(r) = ?- — 1. We denote by F the 
anti-derivative of / which vanishes at 0, and set P{p) = pf{p) — F{p). As pointed 
out before, from a physical viewpoint, P is the pressure. 

The first part of the appendix aims at listing (and deriving formally) the classical 
conservation laws for (jA.2p and (|A.ip . 

Phase invariance. For every a G K, one has 

ip solution of (|A.2p <=> e'"?/) solution of (|A.2p . 

The phase invariance is not seen at the level of (jA.l[) (this amount to changing (j) 
into (j) + £ce)- 

By Noether's theorem or by an easy computation, this leads to the conservation 
of mass: 

(A.3) M := = / pdx. 
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Time translation invariance. For every r € M, one has 

tp{t,x) solution of (|A.2p + t, 2:) solution of (|A.2p . 

By the same time translation, this is expressed at the level of (|A.ip and leads to 
the conservation of the energy (or Hamiltonian): 

(A.4) n:=J (y IV^P + Fi\i,f)) dx = j (ip^l^ + ^p\^ + F{p)) dx. 

Space translation invariance. For every xq E R'^, one has 

^'{t.x) solution to (|A.2p x + xq) solution to (|A.2[) . 

By the same spaee translation, this is expressed at the level of (|A.1[) and leads to 
the conservation of momentum: 

(A. 5) V :^hi\ j etpVilidx — j pvdx. 

Invariance by spatial rotation. Let i? be a spatial rotation. Then 

■il;{t,x) solution of (|A.2|) <^ tpit, Rx) solution oi (IA.2I) . 

By the same spatial rotation on a solution of (|A.1|) this leads to the conservation 
of angular momentum, which we write in the case of R'' for the sake of simplicity: 

(A. 6) A Im J ^ ^ e-ipVij} dx = J x A pv dx. 

Galilean invariance. For every ^0 G 1^'', one has 

■)pit,x) solution of (TOI) e~'^°''°e~'"^l^"l"i/;(t.x + £fnt) solution of (|X2|) . 
For (jA.ip . this implies 

{v,p){t,x) solution ■^=> {v{t,x + e£,ot) — e£,o, p{t,x + e^ot)^ solution, 
and leads to 

dX f f 

(A.7) 'dt^^ with A" := x\tp\'^ dx = xpdx. 

Scale invariance. When /(r) — r" , one has for all A > 0, 

i){t,x) solution of ((X2|) X^l" '\\){}?t,\x) soXwWow oi (|X2|) . 

which implies for (|A.ip 

{v,p){t,x) solution <J=^ (Aw, A^/'^p) (A^i, A.t) solution. 
The associated conservation law is (recall that d stands for the space dimension): 

^ = y"(e2|v*|2 + rfP(|^|2)^dx = 2H + J {dP - 2F){\ij\^) dx 



or, 



dT 
~dt 



j {p\v\^ + dP{p)^ dx = 2H + j {dP - 2F){p) dx 



with 

(A. 8) := Im j e-0x • V?/' = j px-vdx and P{r) rf{r) — F{r). 

Equality (|X8|) remains formally true if the nonlinearity is not a pure power. 
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Momentum of inertia or virial. A direct computation shows that the momentum 
of inertia (or virial) 

(A.9) X:=\j\x\''\i^\^dx=\j\x\^pdx 
satisfies 

dt 

Pseudo-conformal invariance. Let if the pseudo-conformal transform of defined 

by 

e 2et - / e X 

One notices that 



Thus 



For the i^-critical power nonhnearity /(r) = r^/'', one checks that 

-0(i,a;) solution of (|A.2p ■^=^ (^(t, 2;) solution of (|A.2|) . 
For (jA.ip . this yields 



9 / \ 9 

2; V f £ X\ £ \ f £ X 



{v,p){t,x) solution-^ [j ^ \t ) '^vT'T/j solution. 

Using the conservation of energy for ip, one deduces after a lengthy computation 
that 



(A.IO) 



with 



JidP- 2^)(IV'I') dx = 

Z{t) y(i|(x + *£<V)VP + i'i^(|^n)da;, 

= / [\p\^ -tv\' + \/e\^^p\^ + t'Fm')) dx. 

This equality can be proven more simply by using the fact that the operators 

2 

i£dt + and x + i£tW commute. The conservation law for Z is not independent 
from the preceding ones: by expanding the square of the modulus, one observes 
that 

Z{t) =t'^n -tT + 1. 

Consequently, 

dZ ^ dT dl dT 

— ^2tl-L-T - t— + — = 2tl-L- t—, 
dt dt dt dt 

and one recovers (jA.lOp . 
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The Carles- Nakamura conservation law. Formally, the quantity 

U{t) := Re / -0 (x + ieW)^p dx 



is constant. This can be checked by an indirect way (see [17j ) or by a direct 
computation. The interpretation in terms of (jA.l[) is 

d f 
(A.ll) dt^^^ U[t) ^ j p{x-tv)dx. 

This conservation law also results from the conservation of momentum and from 
the conservation law for X since in both cases, (|A.2|) or (jA.l|) . one has 

U^X -tv. 

A. 2. Conservation laws for general capillary fluids. Wc here study to what 
extent the conservation laws listed in the previous subsection are relevant for general 
inviscid capillary fluids. We recall that such fluids are governed by System (j3.7p . 
Throughout, we assume the capillarity k to be a differentiable function on 

For smooth solutions with a non vanishing density. System p.7p recasts in the 
following conservative forn:Q: 



(A.12) 



dt{pv) + div {pv ®v)+ \7P[p) — div K, 
dtp + div(pu) = 0, 



with Kip,yp) (i(K(p) +pK'(p))|Vp2| +pK(p)Ap)/rf - K{p)Vp®Vp. 

Theorem A.l. // (p, w) is a sufficiently smooth solution of (j3.7p which decays at 
infinity, with p non vanishing, then the equalities (|A.3p . (|A.5|) . (|A.6p . (jA.7[) and 
(jA.9|) are still valid, and the energy 

n:= J{lp\v\' + ^\Wp\' + F{p))dx 

is conserved. 

Furthermore, the quantity J- defined in (jA.SP satisfies 
d 



dt 



F=2n + J {dP-2F+ ^(pK)'|V/9p^ dx. 



Proof. It is very likely that most of those conservation laws could be derived from 
Noether theorem. They can also been obtained by pedestrian computations. 

For (jA.3|) . (|A.7p and (|A.9|) . there is no difference with the QHD since the only 
equation on p is concerned. 

For (|A.5|) the proof is obvious from the conservative form (|A.12|) . 

In order to prove the conservation of energy, the simplest method is to take the 

scalar product with pv of the first equation of (|3.7|) . The treatment of the terms 
on the left-hand side of the equation for v is the same as in the compressible Euler 



With the convention (divi^)^ := djK, 
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equation, and one obtains after several integrations by parts and the use of the 
equation for p, 

I / {Ip\v\' + F{p)) = Jpv v(«Ap + i^'lVpp) dx, 

= J dtp(^KAp+^K'\Vp\^^dx, 

= - j K\/p ■ dt'^pdx - ^ J K'c'tplVpp dx, 

To prove (|A.6[) . one writes (using the Einstein summation convention), 

Ai = / Eijkx^pv'' 



with Eijk = if two of the indices are equal, and equal to the signature of (i j k) 
otherwise. 

Using the conservative form of the equation for pv and integrating by parts, one 
thus deduces 

= j e,,x^d,K,,dx- j e.,,x'd,P- j e.,,x^Q,[pv-v^)dx, 

= - j EtjkKjkdx + O- J Eijkpv^v'' dx. 

The tensors K and pv®v being symmetric, the right-hand member of the previous 
equality vanishes. 

The simplest way to prove the conservation law on J- , is to use the second 
equation in (|A.12p . Integrating by parts the capillary term and treating the other 
terms as in the compressible Euler equation one obtains 

= j X ■ div K dx ~ J X ■ (vP + div {pv ® v)^ dx, 

= J x^d,Kijdx + J{dP + p\v\^)dx. 

Using the expression of the tensor K and integrating by parts, one gets 

= -d J (^^{k + pK)\\7pf + pkAp^ dx + J K|Vp|^ 

= ^ l{pnY\Wp\'dx + I K\Wp\^dx, 

from which the claimed equality results. □ 

Remark A. 2. The case of QHD corresponds to (pk)' = 0. For capillary fluids, the 
presence of extra terms in the conservation law for seems relatively harmless if 
p I— pk{p) is an increasing function. One should be able to prove without difficulty 
that with a pressure law such that dP — 2F ^ 0, one has 
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In the case where p i— > is decreasing, we expect the conservation law for T to be 
the key to proving finite time blow-up results similar to those of [42] . 
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